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Abstract. In this work we study a modified Susceptible-Infected-Susceptible (SIS) 

O ' model in which the infection rate A decays exponentially with the number of reinfections 

. , n, saturating after n = I. We find a critical decaying rate ec(/) above which a 

finite fraction of the population becomes permanently infected. From the mean-field 
solution and computer simulations on hypercubic lattices we find evidences that the 
^»-^, upper critical dimension is 6 like in the SIR model, which can be mapped in ordinary 

i-£^ ' percolation. 
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1. Introduction 

Epidemic models have long been studied in the physics community both to infer 
patterns and to develop policies to stop and prevent epidemics on natural environments 
[H El El H] and as a playground to test theoretical ideas [3 El 13 E]. In the simplest 
models, individuals can be either susceptible or infected and, in the latter case, can 
infect others with a given infection rate. If, after some time, the infected individual is 
susceptible to reinfection the model is known susceptible-infected-susceptible (SIS). If, 
otherwise, infected individuals become immune it is called susceptible-infected-recovered 
(SIR) model. Despite the simplicity of these models, they have succesfully modelled a 
varied class of diseases, like dengue [9l [TOl [11], HIV [121 US] and influenza A (HlNl) 

[n [21 [a [15]. 

A fundamental question in the analysis of epidemic models is whether a disease 
will spread through a finite fraction of the population or will it only affect a small 
fraction of the individuals. Also of great importance is the development of efficient 
immunization strategies, both answers depending fundamentally on the underlying 
topology of interconnections between individuals of the population [T6] . 

Although one could create a variety of infection dynamics and network structures, 
an exact mapping of epidemic models onto reaction-diffusion models [6] and subsequent 
development of operator algebras and field-theories [171 UHl HHI show that a small number 
of universality classes might exist, linking different models to the same quantitative 
behavior around the onset of epidemic state [HI [20]. Scaling analysis also reveal that an 
upper critical dimension might exist for each process, above which dimensional effects are 
irrelevant and critical exponents and amplitude ratios are the same as in the mean-field 
approximation (which mimics infinite-dimensional systems, as the number of individuals 
N -^ oo). 

Effects of network topology might also affect the properties of epidemic models: 
while on networks with randomly connected agents there is a finite threshold separating 
epidemic from non-epidemic states, networks with large connectivity fiuctuations, 
such as scale-free networks, lack epidemic thresholds on both SIS and SIR models 
[3 [211 [221 [231 [Ml [251 [26]. 

In this work we study a modified version of the SIS model. In particular, we consider 
a decaying infection rate A, i.e., each individual j that recovers from the infected state 
in a certain time step t decreases his probability to becomes infected again in the next 
time step t -|- 1 in the form Aj(t -|- 1) = eAj(t), where e is a parameter of the model. 
In addition, the decrease of the infection rate occurs only a limited number / of times. 
This type of infectivity has been known to occur both in plants [^ and in the Simian 
Immunodeficiency Virus infection [28]. We study the model on regular lattices and in 
the mean-field regime, and considering the above-mentioned modification, our results 
suggest that the system undergoes a phase transition at critical values e^l) separating a 
phase where the disease reaches a finite fraction of the population (e > ec) from a phase 
where the disease does not spread out (e < ec). We also find that the upper critical 
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dimension is 6, as in the SIR model [HI U\- 

2. Model and Mean-Field Approach 

The population is classified into two classes: Susceptible (S) and Infected (I). The 
transitions between the states S and I occur according to the following automata rules: 

• If individual j is Infected at time t, it is Susceptible at time t + 1 with probability 

a; 

• If, otherwise, individual j is Susceptible at time t, then with probability A it becomes 
Infected at time t + 1 if it is in contact with an Infected individual at the same 
time; 

• Each individual j starts with X{j,t = 0) = Aq, but A depends on how many times 
the individual was infected before, i.e., if an individual j performed the transition 
/ — )■ S in a certain time step t, in the following time step t + 1 the infection rate A 
will be updated in the following form 

A(j,t + l) = eA(j,t) , (1) 

where e is a parameter less than 1 that controls the decaying of the infection rate. 

• In addition, this decay in the infection rate occurs only a maximum number / of 
times, i.e., each individual has a limited capacity to decrease his probability of 
reinfection. 

This dynamics can be relevant to describe diseases such as fiu, for which our chance 
to be reinfected decreases with age. The above rules define a mean-field-like system 
because all individuals interact with all others. Thus, an analytic approach can be 
developed, based on the standard ordinary differential equations for SIS models. We 
can define Sk as the density of Suceptible individuals that have recovered k times from 
the disease, with A^ = e'^Ao being the corresponding infection rate of these individuals. 
Considering that after some time steps (nothing to do with the steady state of the 
system) all individuals will be in two states, I or Si, we can see that the only important 
equation in the evolution of the system is the equation for 5";, i.e., 

§^a./-W. (2) 

Taking into account that in the steady state, i.e., for t — )■ oo, we have that Si = Sstendy 
and dSi/dt = dSstea.dy/dt = 0, Eq. ([2]) give us 

-^steady [Oi ~ € Ao5'stoady j = , (3) 

where we have used the notations S'gteady and /steady to represent the stationary density of 
Susceptible and Infected individuals, respectively, i.e., S'gteady = S(t -^ oo) and /steady = 
I{t -n- oo). There are two solutions of Eq. ([3]): /steady = and /steady = 1 — (a/Ao)e~', 
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Figure 1. Density of Susceptible individuals S{t) (left side, ft-om top to bottom: 
e — 0.80,0.85, ..., 1.00) and density of Infected individuals I{t) (right side, from top 
to bottom: e — 1.00, 0.98, ..., 0.80) as functions of time t for a population of size 
N = 10^, limiting parameter I — 5 and typical values of e. Data are averaged over 200 
realizations. In this case, for e <^ 0.87 the disease does not spread in the system. The 
parameters are: a = 0.05 and Ao = 0.1. 



where we have used the relation S + I = 1. Notice that the nontrivial solution for /steady 
may vanishes for critical values of e given by 

a) ■ (^) 

Using this result, we can rewrite the expression for /steady as 

I 



' steady 



(5) 



In other words, we have a Disease-free phase, where the disease disappears of the system, 
for e < ec(/), with the critical values ec(/) given by Eq. @. On the other hand, for 
e > ec{l) we have an Epidemic phase, where the disease survives and reaches a finite 
fraction of the population. Considering again the relation 5 + / = 1, we can obtain 
a power-law relation between the stationary density of Susceptible individuals and the 
parameters of the system, 

'S'stcady ~ ( T~ ) ^ ' V^) 

After the development of the analytical solution of the problem, we can confront 
it with Monte Carlo simulations. We simulated populations of size A^ = 10^, with 
probabilities a = 0.05 and Aq = 0.1 and different values of e and /. We considered 
that 2% of the individuals are initially Infected in the population, and all results were 
averaged over 200 realizations. Following the rules presented in the beggining of this 
section, the algorithm to simulate the problem is as follows: (i) at each time step, each 
Infected individual j returns to the Susceptible state with probability a; (ii) at the same 
time, each Susceptible individual j becomes Infected with probability Xj if a randomly 
choosen node is Infected. After each transition / — )■ S", the infection rate decreases in 
the form Xj -^ e Xj. 
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Figure 2. Stationary density of Susceptible individuals S'stcady for / = 5 and some 
values of e in the log-log scale. The straight line has slope —5 (left side). For 
the stationary density of Infected individuals /steady we do not have a power-law 
dependency on the parameter e, as predicted analytically in Eq. ([5]). Thus we plot 
here I*^^g^^y = 1 — ^steady, and the slope of the straight line is —5 (right side). The 
parameters are: N = 10^, a — 0.05 and Aq — 0.1. Each point is averaged over 200 
realizations. 
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Figure 3. Density of Susceptible individuals S{t) (left side, from top to bottom: 
e = 0.90,0.92, ..., 1.00) and density of Infected individuals I{t) (right side, from top 
to bottom: e = 1.00, 0.98, ..., 0.90) as functions of time t for N = 10^ and limiting 
parameter I = 10. Data are averaged over 200 realizations. Notice that in this case 
the disease disappears of the system for e <^ 0.93. The parameters are: a = 0.05 and 
Ao = 0.1. 



In Fig. [T]we exhibit results for the density of Susceptible individuals S(t) and the 
density of Infected individuals I{t) as functions of the simulation time t for / = 5 and 
typical values of e. We can observe that the system reaches steady states for all values 
of e, and that for e <~ 0.87 the disease disappears of the system, i.e., we have J = for 
large t. Considering the mean-field calculations derived in the beggining of this section, 
the analytical prediction of Eq. (jl]) for this critical value is e^l = 5) = 0.87055, in 
excellent agreement with the numerical result. 
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Figure 4. Stationary density of Susceptible individuals ^steady for / = 10 and some 
values of e in the log- log scale. The straight line has slope —10. The parameters are 
a — 0.05 and Ao = 0.1. Each point is averaged over 200 realizations. 



We can analyze how the stationary values S'steady = S(t — > oo) depend on the 
parameter e. For this purpose, we have considered time averages of the density of 
susceptible individuals (after the system reach the steady states), and in addition these 
values were averaged over 200 realizations of the system, for each value of e. We can 
observe the power-law behavior (see Fig. |2l left side) 

•S'steady ~ £ '^ , (7) 

with 1/ = 5 in this case (/ = 5). Notice that this result is in agreement with the analytical 
prediction, Eq. IQ. On the other hand, for the stationary values /steady = I(t -^ oo) 
we do not have a power-law behavior, as predicted analytically in Eq. (jS]). However, we 



can analyze the behavior of /, 



steady 



1-/. 



steady ; 



for which Eq. ([5]) give us 



steady 



(8) 



or in other words, L 



steady 



', the same behavior observed for S'steady [see Eq. 



Thus, in Fig. |2] (right side), we exhibit the simulation data for /s^ady versus e for / = 5. 
Fitting data, we obtained /s^ady ~ ^~^ ffl i^ agreement with the analytical result, Eq. 

(ED. 

In Fig. [3] we show results for / = 10 and different values of e. Again, the system 
reaches steady states for all e, but the disease spreads in the system only for e >~ 0.93. 
Observe that this critical value e^l = 10) is greater than the value e^l = 5) ~ 0.87, 
as expected, due to the greater capacity of the individuals to decrease their reinfection 
rates. The analytical result of Eq. (jl]) for this critical value is e^l = 10) = 0.93303, 
again in excellent agreement with the numerical result. 

J Since the behavior of S'steady and ^s*toadv ^^ functions of e is the same, in the following we will analyze 
only the stationary density of susceptible individuals, S'steady 
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Figure 5. Stationary density of Susceptible individuals S'stcady as a function of e for 
different values of I. The straight lines are fittings, which give us S'stcady ~ e^'^'-'-', with 
v{l) = I (left side). It is also shown the exponent i^ as a function of I (right side). The 
squares were estimated from the fittings, whereas the line is the analytical prediction, 
Eq. ^. The parameters in both figures are a = 0.05 and Aq = 0.1. 
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Figure 6. Phase diagram of the model in the plane e versus I, separating the Disease- 
free and the Epidemic phases. The squares were estimated from the simulation data, 
whereas the full (red) line is the mean- field prediction, Eq. Q. The parameters are 
a = 0.05 and Aq = 0.1. 



Considering the stationary values S'stcady for different values of e, we can also observe 
a power-law behavior S'stcady ~ ^~^ (see Fig. H]), but now with a different exponent, 
z/ = 10. In other words, we have the general form 

^steady ~ e'^^ , (9) 

with z/(/) = /, which is supported by numerical results for other values of / (see Fig. |5l 
left side). In addition, the analytical prediction of Eq. (|6]) give us the same behavior of 
the above numerical result, Eq. ([9]), as we can see in Fig. |5] (right side). 
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Figure 7. Density of Infected individuals I{t) as a function of time t for the 
model defined on a square lattice. The parameters are I — 2 (left side, from top 
to bottom: e = 1.00, 0.98, 0.96, ..., 0.86) and I = 5 (right side, from top to bottom: 
e — 1.00,0.99,0.98, ...,0.90) and typical values of e. Data are averaged over 200 
realizations. We have used in these simulations L = 1000, a — 0.05 and Aq = 0.1. 



Taking into account the numerical results for a = 0.05, Aq = 0.1 and different 
values of the parameters e and /, we show in Fig. |6] the phase diagram of the model 
separating the Disease-free and the Epidemic phases. The squares are the critical values 
€^{1) estimated from the simulation data, whereas the curve is the analytical prediction 
of Eq. (jlj). We can observe an excellent agreement between the analytical and the 
Monte Carlo results. 

In the next section we will analyze how the presence of a regular topology affects 
the system and the mean-field results presented in this section. 



3. Simulations on regular d-dimensional lattices 

In this section we will analyze the same model presented in the last section, 
but now it will be defined on regular d-dimensional lattices of linear sizes L. The 
algorithm to simulate the problem is as follows: (i) at each time step, each Infected 
individual j returns to the Susceptible state with probability a; (ii) at the same 
time, each Susceptible individual j becomes Infected with probability m A(j, t)/z, where 
m = 0,1, ...,z is the number of infected nearest neighbors of the individual j, z is the 
coordination number of the d-dimensional lattice [z = 2d) and X{j,t) is the infection 
probability of the individual j at a certain time step t. After each transition / — > 5, the 
infection rate of an individual decreases in the form A — )■ e A. 

Initially, we will consider a square lattice. Thus, we have N = L"^ individuals, and 
each individual interact with four neighbors {z = 4). We have observed that the results 
do not depend strongly on the lattice size. Thus, we simulated populations of size up to 
N = 10^ individuals, i.e., for linear sizes up to L = 1000, with probabilities a = 0.05 and 
Aq = 0.1 and different values of e and /. We have considered that 2% of the individuals 
are initially Infected in the population. All results were averaged over 200 realizations. 
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Figure 8. Stationary density of Susceptible individuals S'stoady as a function of e in 
the log-log scale for the model defined on a square lattice. The parameters are I — 2 
(left side) and I = 5 (right side). Notice that the pure power-law behavior of Ec^. © 
is not observed in the 2D case. However, there are two distinct power-law behaviors 
(full and dashed lines). We have used in these simulations L = 1000, a = 0.05 and 
Ao = 0.1. Each point is averaged over 200 realizations. 
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Figure 9. Stationary density of Susceptible individuals 5'stcady as a function of e in 
the log-log scale for the model defined on a simple cubic lattice. The parameters are 
/ = 2 (left side) and I = 5 (right side). Notice that the pure power-law behavior of 
Eq. ([6|) is not observed in the 3D case. However, there are two distinct power-law 
behaviors (full and dashed lines). We have used in these simulations L = 50, a = 0.05 
and Ao = 0.1. Each point is averaged over 200 realizations. 



In Fig. [7] we show results for the density of Infected individuals as a function of 
time for I = 2 (left side) and / = 5 (right side). We can see that, at least qualitatively, 
the results are the same as in the mean-field limit. However, the critical values edl) 
are different. As examples, we can see in Fig. [7] that we have e^l = 2) ~ 0.88 
and ec(/ = 5) ~ 0.95, whereas we have for the mean-field case e^l = 2) ~ 0.71 and 
ec(/ = 5) ~ 0.87. In other words, the difference bewteen the two analyses increases for 
increasing values of /. This is a consequence of the presence of a topology (neighborhood) 
in the model. 
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Figure 10. Stationary density of Susceptible individuals S'stcady as a function of e in 
the log-log scale for the model defined on hypercubic lattices with dimensions d — A 
(left side) and d = 6 (right side). Notice that the pure power-law behavior of Eq. ([5]) is 
observed for d > 4, and that for d = 6 we have essentially the same behavior as in the 
mean-field case, i.e., we have S'stcady ^ e~'- The parameters used in the simulations 
are a — 0.05, Aq — 0.1, and the lattice sizes considered were L — 20 (for d — 4) and 
L = 8 (for d = 6). Each point is averaged over 200 realizations. 



Following the procedure of the previous section, we can analyze the behavior of 
the stationary density of Susceptible individuals S'steady as a function of e, for different 
values of /. We show in Fig. [8] results for / = 2 (left side) and / = 5 (right side). We can 
observe deviations of the single power-law behavior given by Eq. ([6]). The S'steady values 
follow power laws with two different exponents: one for values of e near the critical point 
(bigger slope, dashed lines in Fig. [8]), and another to intermediary and large values of 
e, with a smaller slope (full lines in Fig. |H]). The same behavior was observed for the 
system defined on a simple cubic lattice {z = 6), as we can see in Fig. |9l This may 
be viewed as a consequence of a small number of neighbors. In fact, if we consider 
higher-dimensional lattices we recover a similar mean-field behavior, i.e., a power-law 
dependency of the S'steady on the parameter e (see Fig. [10]). Considering for example 
/ = 5, the numerical results give us S'steady(^ = 5) ~ e^^-'^'^, S'steady(^ = 5) ~ g-5.55^ 

S'steady (^ = 5) ~ e~^'^'^ and S'steady (^ = 5) ~ g-5.34 ^qj- (^ = /^^ ^ = ^^ (^ = Q g^^ld d = 7, 

respectively. In other words, these results suggest that for d > 6 the system presents a 
similar behavior observed in the mean-field level, with a difference less than 7% to the 
mean-field exponent [see Eq. (J6])]. A similar behavior was also observed for other values 
of the parameter / as / = 2 and / = 10. 

This picture becomes more clear if we estimate the critical values e^l) for different 
dimensions d. In Fig. [TT] we exhibit the phase diagram of the model for some values of 
d. Notice that the Epidemic phase increases for increasing values of the dimensionality 
d. This is a consequence of the increasing number of neighbors (or the coordination 
number z) on higher- dimensional lattices: it is easier to infect an individual that has 
more neighbors. It can be also observed in Fig. [11] that the critical values of edl) for 
d = 6, d = 7 and for the mean-field case are indistinguishable, which reinforces that the 
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Figure 11. Comparative pliase diagram of tlie model in the plane e versus I for 
different dimensionalities d. The symbols were estimated from the simulations. Notice 
that the Epidemic phase increases for increasing values of d, and that the critical 
values ec{l) for d > 6 and for the mean-field approach are the same. The parameters 
are a = 0.05 and Aq = 0.1. 



upper critical dimension of the model is d = 6, as in the SIR model [6l [7]. This is due 
to the frozen state of the individuals with very low infectivity that results from multiple 
reinfections. 



4. Conclusions 

In this work we studied a modified Susceptible-Infected-Susceptible (SIS) model in 
which we have considered that each individual in the population that recovered from 
the disease decreases his probability of reinfection. This decrease occurs a maximum 
number I of times for each individual. This dynamics can be relevant to describe diseases 
such as flu, for which our chance to be reinfected decreases with age. 

Firstly, we have analyzed the problem in the mean-field limit. In this case, 
every individual interact with all others, and we studied the problem with numerical 
simulations and analytical calculations. Considering the initial infection rate Aq (for 
the transition S* — ^ /) and the recovering rate a (for the transition / — )► S*), we have 
found a power-law dependency between the stationary density of susceptible individuals 
5'steady and the parameter e that controls the decaying of the infection rate in the form 
'S'steady = (a/Ao)e~^ In addition, we showed that the system undergoes a phase transition 
at critical values ec(/) = (a/XoY^^ separating a phase where the disease reaches a finite 
fraction of the population (for e > ec) from a phase where the disease does not spread 
out (for e < ec). All results were confirmed by Monte Carlo simulations. 

Considering d- dimensional regular lattices, we have studied the model only with 
numerical simulations. The evolution of the density of Susceptible and Infected 
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individuals is qualitatively similar to the evolution obtained in the mean-field case, but 
the presence of a neighborhood modify some characteristics of the model. In particular, 
the critical values edl) are greater than in the mean-field case for d < 6, which implies 
that the epidemic phase decreases for decreasing values of d. However, for d = 6 and 
d = 7 we have obtained the same values of the critical points edl) as in mean-field 
calculations. The stationary density of susceptible individuals S'steady depends on the 
parameter e in a power-law form only for d > 3, whereas for d = 2 and d = 3we have the 
combination of two power laws. These numerical results suggest that the upper critical 
dimension of the model is rf = 6 as in the SIR model. This is possibly due to the frozen 
state of the individuals with very low infectivity that results from multiple reinfections. 
Thus, the mean-field analytical calculations describe qualitatively well the model, 
providing us the phenomena that are expected to be observed, i.e., it predicts the phase 
transition and the power-law dependency between some quantities of interest. However, 
it fails quantitatively, as it is common in mean-field approximations, because predicts 
different values of the critical points and different power-law exponents. 

Acknowledgment s 

The authors acknowledge financial support from the brazilian funding agency CNPq. 

References 

[1] A. M. Smith, A. S. Perelson, Wiley interdisciplinary Reviews: Systems Biology and Medicine 3, 

429 (2011). 
[2] P. Baccam, C. Beauchemin, C. A. Macken, F. G. Hayden, A. S. Perelson, Journal of Virology 80, 

7590 (2006). 
[3] B. Gonci, V. Nemeth, E. Balogli, B. Szabo, A. Denes, Z. Kornyei, T. Vicsek, PLoS ONE 5(12): 

el5571 (2010). 
[4] R. M. Anderson, R. M. May, Infectious Diseases of Humans: Dynamics and Control (Oxford 

University Press, 1991). 
[5] J. Marro, R. Dickman, Nonequilihrium Phase Transitions in Lattice Models (Cambridge University 

Press, Cambridge, 1999). 
[6] H. Hinrichsen, Advances in Physics 49, 815 (2000). 
[7] M.E.J. Newman, Phys. Rev. E 66, 016128 (2002). 
[8] N. T. J. Bailey, Ttie Mathematical Theory of Infectious Diseases and its Applications (Hafner 

Press, New York, 1975). 
[9] M. Derouich, A. Boutayeb, E. H. Twizell, BioMedical Engineering 2:4 (2003). 
[10] L. Esteva, C. Vargas, J. Math. Biol. 46, 31 (2003). 

[11] N. Nuraini, E. Soewono, K. A. Sidarto, Bull. Malays. Math. Sci. Soc 30, 143 (2007). 
[12] R.M.Z. dos Santos, S. Coutinho, Phys. Rev. Lett. 87, 168102 (2001). 
[13] F. Baryarama, L. S. Luboobi, J. Y. T. Mugisha, Amer. J. Infect. Diseases 1, 55 (2005). 
[14] I. Gordo, M. G. M. Gomes, D. G. Reis, P. R. A. Campos, PLos ONE 4(3):e4876 (2009). 
[15] G. Katriel, L. Stone, PLoS Curr. 1: RRN10460 (2009). 

[16] R. Cohen, S. Havhn, D. ben-Avraham, Physical Review Letters 91, 247901 (2003). 
[17] M. Doi, J. Phys. A 9, 1479 (1976). 
[18] H. K. Janssen, Z. Phys. B 42, 151 (1981). 
[19] P. Grassberger, A. de la Torre, Ann. Phys. 122, 373 (1979). 



Critical behavior of the SIS epidemic model with time- dependent infection rate 13 

[20] P. Grassberger, Z. Phys. B 47, 365 (1982). 

[21] R. M. May, A. L. Lloyd, Phys. Rev. E 64, 066112 (2001). 

[22] Y. Moreno, R. Pastor-Satorras, A. Vespignani, Eur. Phys. J. B 26, 521 (2002). 

[23] R. Pastor-Satorras, A. Vespignani, Phys. Rev. E 63, 066117 (2001). 

[24] R. Pastor-Satorras, A. Vespignani, Phys. Rev. Lett. 86, 3200 (2001). 

[25] R. Parshani, S. Carmi, S. Havhn, Phys. Rev. Lett 104, 258701 (2010). 

[26] M. Boguiia, R. Pastor-Satorras, A. Vespignani, in Statistical Mechanics of Complex Networks, 

edited by R. Pastor-Satorras, M. Rubi and A. Diaz-Guilera (Springer Verlag, Berhb, 2003), vol. 

625 of Lecture Notes in Physics. 
[27] C. A. Gilligan, S. Gubbins, S. A. Simons, Phyl. Trans. R. Soc. Lond. B 352, 353 (1997). 
[28] N. K. Vaidya, R. N. Ribeiro, C. J. Miller, A. S. Perelson, Journal of Virology 84, 4302 (2010). 



